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Abstract. We study the random graph G„_A/n conditioned on the 
event that all vertex degrees lie in some given subset S of the non- 
negative integers. Subject to a certain hypothesis on S, the empirical 
distribution of the vertex degrees is asymptotically Poisson with some 
parameter ^ given as the root of a certain 'characteristic equation' of S 
that maximises a certain function tps{f^)- Subject to a hypothesis on S, 
we obtain a partial description of the structure of such a random graph, 
including a condition for the existence (or not) of a giant component. 
The requisite hypothesis is in many cases benign, and applications are 
presented to a number of choices for the set S including the sets of (re- 
spectively) even and odd numbers. The random even graph is related 
to the random-cluster model on the complete graph Kn- 



1. Introduction 

Let 5 be a fixed nonempty set of non-negative integers. The purpose of 
this paper is to study the structure of random graphs having ah their vertex 
degrees restricted to the set S. 

We call a graph an S-graph if all its vertex degrees belong to S. For 
example, if 5 = {s} is a singleton, an 5-graph is the same as a regular 
graph of degree s. (We are not going to say anything new about this case.) 
One of our main examples is the class of Eulerian graphs, or even graphs, 
given by the set of even numbers S = 2Z>q, with Z>o the set {0, 1,2,...} 
of non-negative integers. See Section [6] for further examples. 

More precisely, we will study the random graph Gn,p-s defined as Gn,p 
conditioned on being an S-graph, where Gn,p is the standard random sub- 
graph of the (labelled) complete graph Kn where two vertices are joined 
by an edge with probability p G (0, 1), and these (2) events, corresponding 
to the edges of K^, are independent. In other words, Gn,p-s is a random 
5-subgraph of Kn such that, if G is any given subgraph of Kn that is an 
5-graph, then 

„c(G)/i_ )(^)-c(G) 

P(G„,p;5 =G) = Z \ P> -, (1.1) 

r{Gn,p IS an 5-graph) 
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where e(G) is the number of edges of G. We are interested in asymptotics 
as n — > oo, and we wih tacitly consider only n such that there exists an 
5-graph with n vertices, in other words, such that the denominator in (jl.ip 
is non-zero. Thus, a finite number of small n may be excluded; moreover, if 
S contains only odd integers, then n has to be even. (It is easy to see that, 
apart from this parity restriction, all large n are allowed.) 

Remark 1.1. The choice p = \ gives a random 5-graph that is uniformly 
distributed over all 5-graphs on n labelled vertices. However, we will in this 
paper instead study the case when p is of order 1 /n and the average vertex 
degree is bounded (for Gn,pi and as we shall see later, for Gn,p-s too). 

It follows immediately from (jl.ll) that two 5-subgraphs of Kn with the 
same degree sequence are attained with the same probability. Hence, the 
conditional distribution of Gn,p-s given the degree sequence is uniform. We 
will therefore focus on studying the random degree sequence of Gn,p-s] it is 
then possible to obtain further results on the structure of Gn,p-s by applying 
standard results on random graphs with given degree sequences to Gn,p-s 
conditioned on the degree sequence. For example, using the results by Mol- 



loy and Reed IJ, ll5| we obtain Theorem 13.11 below on existence of a giant 



component in Gn,p-s- 

2. Main theorem 

By symmetry, the labelling of the vertices and thus the order of the de- 
gree sequence is not important, and we shall therefore study the numbers 
of vertices with given degrees, rather than the degree sequence itself. We 
introduce some notation. 

Let M be the set of sequences n = (no,ni, . . . ) of non-negative integers 
rij £ Z>o with only a finite number of non-zero terms nj. Let 

Afs := {n e Af : nj = when j ^ S}, 

the set of such sequences supported on S. For a (multi)graph G, let nj{G) be 
the number of vertices of degree j in G, j G Z>o, and let n(G) := {nj{G))°^Q 
be the sequence of degree counts. Thus, G is an 5-graph if and only if 
n(G) G Ms- Clearly, cf. (fTT]) . 

P(G„,p is an 5-graph) = /^^^l - p)^^)-^^^). (2.1) 

We shall call this summation the partition function, denoted as Z^.p-^s- 

Let V be the set of probability distributions on Z>o. In other words, V 
is the set of sequences tt = (7ro,7ri, . . . ) of non- negative real numbers such 
that TTj = 1. We regard V as a topological space with the usual topology 

of weak convergence (denoted — ^); it is well known that this topology on 
V may be metrised by the total variation distance 



dTv(7r,7r') := 



/ 1 

7Tj\ 
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If G has n = '^jnj{G) vertices, let T^j{G) := nj{G)/n, the proportion of 
vertices of degree j, and 7r(G) := n{G)/n = (vrj(G))^Q. Note that 7v{G) is 
the probabihty distribution of the degree of a randomly chosen vertex in G. 
Let (j)s be the exponential generating function of S, 



^s(^')■■=Y.^■ (2.2) 

Note that this is an entire function of /x, and that (1)3 {/J-) > for /j, > while 
05(0) > if and only if G 5. 

Let Pos{fi) be the distribution of a Po(/i) distributed variable given that 
it belongs to 5, i.e., C{X^ \ £ S) with ~ Po(^). Thus, recalling 

Posif^m = kGS. (2.3) 

This conditional distribution is always defined for /i > 0, and in the case 
G 5 for fj, = too (in which case it is a point mass at 0). The mean of the 
Fos{n) distribution is 

E(X,|X„eS) = -l-E^ = ^^ (2.4) 



Let A > 0. We shall refer to the equation 



(2.5) 



05 (m) a ' 

as the characteristic equation of the set S (for this value of A), and we write 

where we allow = only if € 5. We further define the auxiliary function 

V'5(^):=log05(^)-f|^ (2.7) 
and note that when /i G E{X), ipsifJ-) equals the simpler function 

V^5,i(/i;A) :=log05(Ai)-|^. (2.8) 

All logarithms in this paper are natural. We let ci, Ci, . . . denote positive 
constants, generally depending on S and A (or (A^)) and sometimes on other 
parameters too (but not on n), which may be indicated by arguments. We 
sometimes assume that n > 1 to avoid trivialities. 

Theorem 2.1. Let A„ ^ A > and suppose that E{X) contains a unique 
fl = /1(A) that maximizes il^s[f^) (or, equivalently, 'ips^i{^; X)) over E{X). 
Then, the following hold, as n —> oo : 
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(i) '^{Gn,\r./n;s) ^ Po5(/I). In other words, for every j G S, 

^P°^(^)W = ^- (2.9) 

(ii) All moments of the random distribution T^{Gn^\^/n;s) converge to the 
corresponding moments of P05 (/I) . In other words, if di, . . . ,dn is 
the degree sequence of Gn^\„/n;S' -^t^ ~ Po(/u), then for every 
re {0,00), 

n,X„/n;S) 



n ^-^ n 



i=l j=0 



(2.10) 



3=0 

In particular, 

._E(X^|X^G5)-^^. (2.11) 

(iii) The error probabilities in (i) decay exponentially: for every e > 0, 
there exists a constant c\ = ci{e,\,S) > such that, for all large n, 

P(dTv(vr(G„,;,„/„.5),Po5(/I)) >e)< e-^^". (2.12) 

(iv) We have that 

^ log^„,A„/n;5 = ^ logP(G„^A„/n IS an 5-graph) il^sCfi) - ^A. (2.13) 

More generally, let Eq{X) be the subset of E{X) where ipsifJ-) (or V'5,i(/^; A)) 
is maximal: 

Eo{X) := { fx e E{X) : ^l;s{Ii) = max V'-S (/"')) • (2.14) 

If Eq{X) contains a single element, we thus take that element as /i(A); in 
particular, if \E{X)\ = 1, then E{X) = Eo{X) = {/i(A)}. We shah see in 
Section m that this is the normal case: E{X) is always finite and non-empty, 
and |£'o(A)| = 1 except for at most a countable number of values of A. 
Further results on fi and the auxiliary functions are given in Sectional 

Remark 2.2. The set Eo{X) may contain more than one element, see Exam- 
ple [HiHl In this case, the theorem may not be applied, but the proof in Sec- 
tions EHS] extends to show that the random degree distribution Tr{Gn^\^/n;s) 
approaches the finite set F := {Po5(/u) : /j, G -E'o(A)} in the sense that the 
analogue of (j2.12p holds for the distance to this set, i.e., 

P( min dTv(7r(G„,;,„/„.5),Po5(/i)) > e) < fi-^^". (2.15) 

We can regard the distributions in F as pure phases, in analogy with the 
situation for many infinite systems of interest in statistical physics, but for 
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the finite systems considered here this has to be interpreted asymptotically. 
Thus, for large n, the degree distribution of Gn,\„/n;S is approximately given 
by one of the pure phases, but we do not know which one. It follows that if we 
let n oo, one of the following happens for the random degree distribution 
TT = Tr{Gn,\„/n;s) (regarded as an element of V): 

(i) vr converges in probability to Po5(/i) for some fj, G Eq{X). 

(ii) vr converges in distribution to some non-degenerate distribution on 
F. (A mixture of two or more pure phases.) 

(iii) There are oscillations and vr does not converge in distribution; suit- 
able subsequences converge as in (i) or (ii), but different subse- 
quences may have different limits. 

It is easy to show by a continuity argument that all three cases may occur 
in Example 16.81 for suitable sequences (A„) (with A„ Xq defined there). 
We do not know whether all three cases may occur for fixed A. 
We shall not investigate the case |£'o(A)| > 1 further here. 

Remark 2.3. The reason for taking n large in (j2.12p is as follows. Suppose, 
for example, that S is an infinite set. Then the left hand side of (j2.12p 
trivially equals 1 for any fixed n and sufficiently small e. One way around 
this, at least when A„, = A, is to replace the right hand side of (I2.12p by 
2e-<=2n foj, suitable C2 > 0. 

We close this section with an informal explanation of the results of Theo- 
rem I2.lt several applications of which are presented in Section [6l Recall the 
partition function Z^^x/n-s of (|2.ip . considered as a summation over suit- 
able graphs. We wish to establish which graphs are dominant in this sum- 
mation. In so doing, we will treat certain discrete variables as continuous, 
and shall study maxima by differentiation and Lagrange multipliers. Let 
TT = (vro, vri, . . . ) be a sequence of non-negative reals satisfying vrj = 1, 
and vrj = for i ^ S. We write 

h' = i^(7r) = ivrj. 

i 

Let Z(7v) represent the contribution to the summation of (12.11) from graphs 
G having, for each i, approximately nvrj vertices with degree i. The (empir- 
ical) mean vertex-degree of such a graph is ly. 

Now, Z{7r) is a summation over simple graphs subject to constraints on 
the vertex degrees. It may be approximated by a similar summation Z'{7v) 
over certain multigraphs, and this is easier to express in closed form, as 
follows. The number of ways of partitioning n vertices into sets Vo,Vi, . . . 
of respective sizes nvr^, i > 0, is 

n! 

(nvro)!(nvri)! • • • 

Each vertex v £ Vi will be taken to have degree i, and we therefore provide v 
with i 'half-edges'. Each such half-edge will be connected to some other half- 
edge to make a whole edge. Since half-edges are considered indistinguishable. 
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we shall require the multiplicative factor 

The total number of half-edges is 2N = n^ - ivTj = nz^, and we assume for 
simplicity that A*" is an integer. These half-edges may be paired together in 
any of 

(2iV-l)!! = (2iV-l)(2iV-3)---3-l = ^ 
ways, and each such pairing contributes 



to Z'(7v). We combine the above to obtain an approximation to Z{tt): 

By Stirling's formula, as n ^ oo, 

-logZ(7r) ^ \v\og{\v/e) - - ^ ^aog(i! vTi). (2.16) 

We maximize the last expression subject to vTj = 1 to find that 

^^ = A^, iGS, (2.17) 

for some constant A and some fj, satisfying 

fi = (2.18) 

Thus TT is the mass function of the Po5(//) distribution and, by ()2.4p and 
the definition of 



(2.19) 



We combine this with (|2.18|) to obtain the 'characteristic equation' 

If there exists a unique fi satisfying the characteristic equation, then we 
are done. If there is more than one, we pick the value that maximizes the 
right hand side of (|2.16|) . That is to say, the exponential asymptotics of 
^n,x/n;S ^^'^ dominated by the contributions from graphs with tt satisfying 
(j2.17p with /i chosen to satisfy the characteristic equation and to maximize 

Note from (pT6]) that 

-logZ(7r)^V5(Ai)-^A, (2.20) 

n 

and part (iv) of Theorem (j2.ip is explained. 
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We make the above argument rigorous in the forthcoming proof of Sections 
[IHHl a substantial part of which is devoted to proving that the conditional 
distribution of vertex degrees is concentrated near its mode. 

3. The giant and the core 

We show next how to apply Theorem 12. 1^ in co njun ction with results of 
Molloy and Reed 14, [l^ and Janson and Luczak [lH, to identify the 



sizes of the giant cluster and the /c-core of Gn^\^/n-s- The proofs are deferred 
to Section [TOl 

We consider first the existence or not of a giant component in the random 
5-graph Gn^x„/n;S as A„ ^ A > 0. Let tt = (ttq, tti, . . . ) be a vector of non- 
negative reals with sum 1, and write = Yljj'^j- explained in (l^: [l5|. 
if we consider the random graph with given degree sequence d = (c?j)", and 
assume that there are n(l + o(l))7rj vertices with degree j, the quantity that 
is key to the existence of a giant component is 

j 

Subject to certain conditions, if Q^tv) > 0, there exists a giant component, 
while there is no giant component when Q{7v) < 0. 

We shall apply this with ttj = Pos{li){j} , and to that end we introduce 
some further notation. Let, see (j2.4p . 



uif.):=Y.jFos{f.){j} = i^, (3.1) 



Qif^) ■■= Y.^{j - 2) po.(^){j} = (3.2) 



Note that XAt(O) = X^iX) = 0- Furthermore, the only possibly negative term 
in the sums in (j3.2p and (|3.3p for ^ G (0, 1) are those with j = 1, while the 
terms with j = and j = 2 always vanish and the others are positive unless 
Po5(/x) = 0. 
Let 

Fn,p;S be the component of Gn,p;S with the largest number of vertices, 
and let L^ ^.^ be the second largest. (Break ties by any rule.) We write v(i7) 
for the number of vertices in a graph H. 

Theorem 3.1. Suppose that S ^ {0,2}. Let A„, ^ A > and suppose 
that Eq{\) contains a unique element > 0. Then, G^ Xn/niS ^ giant 
component if and only if Q{fl) > 0, i.e., if and only if ^4)g{'jX) > 4>'g{Ji). 
More precisely, as n ^ oo, 



n 
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where 



05 (^) 

C = C(/i):=H/^)(l-?) (3.5) 

with ^ = G [0, 1] given as follows: 

(i) > and 1 G 5, i/ien ^ G (0, 1) is i/ie unique solution to 
X-fiiO = with < ^< 1, and 7, C > 0; 

(ii) if Q{fi) > and I ^ S, then ^ = and 7 = 1 - 1/0^ (/I) > 0, 
C = ^i^iJl) > 0; ^ 

(iii) if Q{'iT) < 0, i/ien ^ = 1 and 7 = ^ = 0. 

Remark 3.2. If /2 = 0, then = and we are in Case (iii) with no 

giant component; in fact, by Theorem l2.1l uq/u — ^ 1 so ahnost all vertices 
are isolated. In this case Xjii^) = for all ^. 

If /I > 0, then < on (0, ^) and x/i > on 1) in all three cases, as 
follows from [Ilj, Lemma 5.5], which yields another characterization of ^. 

Remark 3.3. If 1 ^ 5, then Q{'jl) > as soon as /i > 0. Hence we are in 
Case (ii) if /I > and in Case (iii) in /I = 0. 

Remark 3.4. We have excluded the cases 5 C {0,2}, i.e., the trivial case 
S = {0} and the cases {2} and {0, 2} that are exceptional; in the latter 
cases n~^\{Tnx^^in.g) has a continuous limiting distribution and thus not a 
constant limit when /i > 0, see Examples 16.21 and 16.71 We note also that 
Q(^) = for all ^ in the excluded cases. 



Remark 3.5. It is easily seen, using (13. 3p . that ^ equals the extinction 
probability of a Galton- Watson process with offspring distribution 

jP05(/I){i} 



\X = j - 1, 

V 



iG5\{0}, 



that is, the distribution Po5_i(/i) where S — 1 := {k > : k + 1 G S}. (Note 
that (/)5_i(//) = Hence 7, the asymptotic relative size of r„,A„/n;5) 

equals by (j3.4p the survival probability of a Galton-Watson process with 
offspring distribution Vos~ 1 (/i) and initial distribution P05 (/i) . 

The /c-core of a graph G is the largest induced subgraph having mimimum 
vertex degree at least k. The A:-core of an Erdos-Rehyi random graph has 
attracted much attention; see [13] and the references therein. Theorem 1 2.11 
may be applied in conjunction with Theorem 2.4 of Janson and Luczak |10l | 
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(k) 

to obtain the asymptotics of the k-coie of Gn,\„/n\S- Let K^^.^ denote the 
fe-core of Gn,p;S- We shall require some further notation in order to state 
our results for K^'^l , 

Let k £ {2,3, . . .}. Let /U > 0, and let be a random variable with the 
Fos{iJ,) distribution. For r G [0,1], let VF^,r be obtained by 'thinning' 
at rate 1 — r so that, conditional on W^, Wfj,^r has the binomial distribution 
Bin(W^,r). For A; G {2,3, . . . }, let 

oo 
l=k 

Theorem 3.6. Let — > A > and suppose that Eq{X) contains a unique 
element 'jl. Let k >2, and let, with v = as above. 



sup{r < 1 : z^r^ = /i^^fc(r)}. 



As n ^ oo: 
(i) if?=0. 



\ n,\„/n;S' ' „ n,\„/n;S> 



if, further, k > 3, then 



(ii) ifr> 0, and in addition vr"^ < h'p_^}^{r) on some non-empty interval 
(f — e, r), then 



1 



n 
n 



e{K [ I c) — > —hrrhir) = —ur' 

^ n,X„/n;S' 2 "'^ ' 2 



Remark 3.7. Let be the Galton- Watson process with offspring distribu- 
tion Po5_i(/x), started with a single individual o, and let be the modified 
process where the first generation has distribution Po5(;u), cf. Remark 13.51 
It may be seen that r is the probability that the family tree of Xfj_ contains an 
infinite subtree with root o and every node having k — 1 children. Similarly, 
h-fi^k{r) equals the probability that X ^ contains an infinite fc-regular subtree 
with root o (the root has k children and all other vertices have A; — 1). It is 
easy to see heuristically that this yields the asymptotic probability that a 
random vertex belongs to the /c-core, see Pittel, Spencer and Wormald [3] 
(for S = Z>o), but it is difficult to make a proof based on branching process 



theory; see Riordan [17| where this is done rigorously for another random 
graph model. 
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4. Roots of the characteristic equation 

To avoid some trivial complications, we assume throughout this section 
that S ^ {0}, thus excluding the trivial case S = {0} for which Gn,x„/n\S 
comprises isolated vertices only. 

Lemma 4.1. < Ci(j)s{fJ,) for all fj, > 1. 

Proof. By 

We split the sum into two parts. For j < 4/i, 



For j > 4^, Stirling's formula implies 



and thus, for /i > 1, 



. oo 

E ^<E4(i)'=^2 = C305(l)<C305(/^). □ 
ie5,i>4Ai j=l 

Theorem 4.2. For each A > 0, the set E[\) is finite and non-empty. 

Proof. The characteristic equation (j2.5p may be written as /i(^) = where 
h{fi) = Xfj.cj}'g{fi) — fi'^cpsifJ')- By Lemma I^Tl for fi > I, 

h{fi)<{CiXfi-fi^)4>s{fi), (4.1) 

and thus /i(/u) < for fi > C4 = maxjCiA, 1}. 

Since h is an entire function, and does not vanish identically by what we 
just have shown, it has only finitely many zeros in each bounded subset of 
the complex plane, and in particular in the interval [0, C4]. Hence E{\) is 
finite. 

To see that E(\) is non-empty, let s be the smallest element of S. If 
s = 0, then G E{X). If s > 0, then h{n) ~ Xsj^^/sl as /i ^ 0, so /i(//) > 
for small positive fi. Since further h{fi) is negative for large fi, h possesses 
a zero on the positive real axis. □ 

We have defined Ji as the maximum point of tps or ips,! on E{X). The 
next theorem shows that, alternatively, it can be defined as the maximum 
point of ^^5,1 on [0,00) (but not of ips)- Furthermore, instead of V's.ij we 
can use the function 

that arises as follows. In Section [71 we will indicate the use of multigraphs 
in proving Theorem 12. 1^ of which we shall derive a multigraph equivalent 
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at Theorem 17.31 With ^^i^.g denoting the multigraph partition func- 
tion, we shah see in the proof of Theorem 17.31 that V'cS,2(/^; A) represents the 
contribution to log Z* ^^^.^ from multigraphs with degree distribution 
close to Po5(/i), see Remark 18.31 For this reason, ^5,2 is a more natural 
function than V'.s.i; although it has a more complicated formula. We shall 
have to exclude the trivial case when S = {s} is a singleton; in this case 
V'<S,2(a*; A) = ^s(log(As) — l) — log(s!) is constant. 

It is easily seen that ips,2{fJ'] ^) ^ V'5,i(A')A) for all fi and A > 0, with 
equality if and only if £ E{X). 

We regard V^.i ^-nd iIjs,2 as functions of /i, with A considered a fixed 
parameter. These functions are evidently analytic on (0,cxd). Note that if 
G 5, then (psi^) = 1 and il^si^t)-, V'<s,i(a';A) and V'-S,2(a';A) are continuous 
at // = with "05(0) = V'5,i(0;A) = "05,2(0; A) = 0. On the other hand, 
if ^ 5, then 05 (0) = and ■05,i(/u; A) —oo while a simple calculation 
yields '05,2(Ai; A) — > ^s(log(As) — l) — log(s!) as ^ — > 0, where s = min5. 

Theorem 4.3. The following hold for every fixed A > and j = 1 or 2, 

except for j = 2 in the trivial case \S\ = 1. 

(i) E{\) is the set of stationary points ofipsj, possibly with added: 

E{X) n (0, oo) = j/^ : ^^sAi^-^ A) = 

(ii) Eq{X) is the set of global maximum points ofipsj- 



Eo{X) = <fi: ilJs^jilJ'] A) = niaxV'5,i(/^'; A) 
L /i'>o 

Proof, (i): Differentiation yields 

and, after some simplifications, 

^^.,2(/.;A) = -(^jlog(-^j. (4.4) 
By the Cauchy-Schwarz inequality, provided /i > and \S\ > 2, 

dii V 05 (/^) ) dix\ J2j(zs l^/j- 



, > 0. (4.5) 

{E,,si^^m" 

(See Theorem 15 . 21 below for a more general result.) By (j4.3p - ()4.5p . for /i > 0, 
ip'gj{fi;\) = if and only if (j)'g{fi)/(j)s{fi) = /u/A, i.e., the characteristic 
equation (12. Sh holds. 

(ii): By (|4.3l) " (j4.5l) and Lemma [AA\ ipsj is decreasing for large /i. Fur- 
thermore, by the remarks prior to the theorem, tpsj is either continuous 
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at or tends to — oo there. This imphes that ip^j has a finite maximum, 
attained at one or several points in [0, oo). It remains to show that the 
maximum points belong to E{X); it then follows that E(){X) equals the set 
of maximum points. 

If > is a maximum point of 4's then ip'^ = and fi £ E{X) by 

(i). 

If is a maximum point and G 5, then E -E'(A) by definition. Finally, 
if ^ 5, and s := min5 > 0, then 4>'g{fi)/(f>s{n) >ooas/x— >0, 
and it follows from ()4.3p ~ (|4.5p that ^jJg j > for small hence is not a 
maximum point in this case. (In fact, for j = 1, "05,1(0) = — oo when 0^5, 
as remarked above.) □ 

We define /i*(A) := min£'o(A) and /i*(A) := max£'o(A); thus ii?o(A)| = 1 
(and Theorem 12.11 applies) if and only if fl^ = /i*, and in that case = fl^ = 
fl*. We have defined /2 only when |£^o(A)| = 1; for convenience we extend the 
definition to all A > by letting /1(A) by any element of -E'o(A) (for example 
/I* (A) or /I*(A)). 

Corollary 4.4. For every A > and j = 1, 2, 

05(/^(A)) = ipsjCfiW; A) = maxV'5,j(/^; A). 

Theorem 4.5. // < Ai < A2, then /2(Ai) < /i(A2), with equality only if 
/J(Ai) = /i(A2) = 0. 

Proof. If /i < /I*(Ai) G So(Ai), then V5,i(/^;Ai) < ^5,i(/J*(Ai); Ai) by The- 
orem [4]3]^ii) , and thus 

'^Ps,l{^J'■A2) = V'5,i(/^; Ai) + Yijl~J^ 

< V'5,i(/S*(Ai); Ai) + ( 1 - 1) = ^5,i(/i*(Ai); A2), 

so is not a global maximum point of 05,i(/u; A2) and, by Theorem 14. 3r iil 
again, fj. ^ £'o(A2). Hence, /i(A2) > /i*(A2) > /i*(Ai) > /u(Ai). Equality 
is possible only if fi := fi{Xi) = /J(A2) S -E'(Ai) R £'(A2), and then the 
characteristic equation (12. 5p is satisfied with and both Ai and A2; hence 
n'^/Xi = /-i^/A2 and fi = 0. □ 

Theorem 4.6. (i) For every A > 0, /i(A') / /I* (A) as A' / A and fi{X') \ 
/I* (A) as X' \ X. 

(ii) /2(A) ^ as A ^ 0. 

(iii) /x(A) — > 00 as A ^ 00. 

Proof, (i): Let /^o := lin^A' /'A/2(A'); the limit exist by the monotonicity in 
Theorem 14.51 For any fixed //, '05,i(/i;A') < ^5,i(/I(A'); A') by Theorem 14.31 
and it follows by continuity that V'5,i(/^; A) < '05,i(/io; A). Hence, by Theo- 
rem 1131 again, /io G -E'o(A), and Theorem 14.51 implies that /io = /U=k(A). 
The second statement is proved similarly. 
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(ii): This is similar. Assume /Uq := limx^Qjl{X) > 0, and let fii := 
Then ijjs,i{fJ'i', A) < tps,i(fi{\); A) for all A by Theorem 14. 31 which contradicts 
the fact that, by 



V'^.ilw; A) - V5,i(m(A); A) = ^^^^^3^ + 0(1) 

>^^ + 0(l)-oo 

as A ^ 0. 

(iii): Assume fiQ := limA^oo m(A) < oo, and let fii := fiQ + 1. Then 
'^s,i{fJ-i'i A) < V'5,i(/^(A); A) for all A by Theorem l4.3l and it follows from (|2.8p 
by continuity, letting A — > oo, that log4>s{ni) < log(f)s{fiQ), a contradiction 
since (ps is strictly increasing. □ 

Remark 4.7. For large A, we have the estimates csX^^"^ < /1(A) < C5A, 
where the lower bound follows from (|2.5p and the upper from (j4.ip . Exam- 
ples 16.11 and 16.61 show that both these orders of growth can be attained. 

We see from Theorem 14.61 that jE'olA)! > 1 exactly when /u(A) is discon- 
tinuous, and that all discontinuities are jump discontinuities: /i jumps form 
/i*(A) to /i*(A). In accordance with Remark 12.21 we interpret these discon- 
tinuities as phase transitions of Gn^\„/n;S- More generally, we say that we 
have a phase transition at each A where ju is not analytic. (See, further. 
Theorem 14.151 ) We show that there is only a countable number of phase 
transitions with jump discontinuities, and we note from Example 16.91 that 
the number may be infinite. 

Theorem 4.8. The set A := {A > : |£'o(A)| > 1} = {A : /I*(A) < /I*(A)} 
of A such that Theorem \2.1\ does not apply is at most countable. 

Proof. By Theorem 14.51 the open intervals (/i=K(A), /i*(A)), A > 0, are dis- 
joint, and thus at most a countable number of them are non-empty. □ 

It follows from Theorem 14.51 and its corollaries that is the inverse 
function of a continuous non-decreasing function with graph {(/i,A) : yu G 
[/2*(A), /2*(A)]}; the exceptional set A consists of the values taken by this 
function on the intervals where it is constant. 

For /i > we can rewrite the characteristic equation (|2.5p as 



A = A(,.):=i^M^>. (4.6) 

Thus, ^(A) = {/X > : A(/i) = A} or {/i > : X{fx) = A} U {0}. Note 
that our assumption S ^ {0} implies that (p'^ifi) > for all /z > 0, so A is 
well-defined. 
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Lemma 4.9. The function A is positive and analytic on (0, oo), with lim^^oo X{f^) 
oo and 



lim A(/i) 



0, 0^5, 

0, Og5, iGcS, 

1, G 5, 1 ^ 5, 2 e cS, 
oo, G 5, 1 ^ 5, 2 ^ 5. 



Proof. That A is analytic and positive is evident. By Lemma [4 .11 X{^) > C4/i 
for large //, and the behaviour as /i — > follows by looking at the first non- 
zero terms in the Taylor expansions of i;^^ and (/>^. □ 

Lemma 4.10. ip'^i/J-) and A'(^) have the same sign for every /x > 0. Hence 
tps and A have the same stationary points and are increasing or decreasing 
on the same intervals. 

Proof. By (|2.7p and (|2.8p . = i^s,li^J''■, X{^)). Differentiating and using 

(j4.3p we obtain, for all /i > 0, 

o o 2 

i^M = ^V'^.il/x, A(/x)) + A(/.))A'(^) = + □ 

Theorem 4.11. (i) //A is decreasing on an interval (/ii,/i2) with < /ii < 
//2, i/ien i/iere exists A G A with /i*(A) < //i < //2 < 

(ii) A = i/ and on/|/ i/ A is increasing on (0, oo). In this case, for 
every A > either E{X) = {/2(A)} with /1(A) > or £;(A) = {0,/2(A)} wii/i 
/2(A) > 0. 

Proof, (i): Suppose that /1*(A) G (/ii,/i2) for some A > 0. Taking a se- 
quence A„ ^ A, we have fl{Xn) /2*(A) by Theorem 14.61 so, for large n, 
/ii < /2(A„) < /i=„(A) < [i2 and hence A„ = A(/i(A„)) > A(/2*(A)) = A, a 
contradiction. Hence, fi^{X) ^ {fii,fi2) for all A > 0. 

Let Ao := sup{A : /I*(A) < /xi}. Thus, /2*(A) > fi2 for A > Aq. By 
Theorem I4.6r ii)(iii). < Aq < oo, and by Theorem 14.6^ 1) /U*(Ao) < fJ-i and 
/2*(Ao) > /t2- In particular, /I*(Ao) < /i*(Ao) so Aq G A. 

(ii): If A is not increasing, then A'(/i) < for some /i > and (i) applies 
to some interval (fj, — e, fi + e) and shows that A 7^ 0. 

If A is increasing, then ips is (strictly) increasing on (0, 00) by Lemma [4.10t 
if G 5, then -0^ is continuous at and thus increasing on [0, 00) also. 
Consequently, -E'(A) contains a unique fi = maxE{X) that maximizes ips. 
Further, when A is increasing, there is at most one positive solution to 
A = A(/i), and thus to (j2.5p . and the result on -E'(A) follows. □ 

We next study whether ju = is possible. Note that this is a rather 
degenerate case, when Theorem 12. II shows that G„,A„/n;5 is very sparse with 

Op(n) edges and uq/h 1, which is to say that n(l — Op(l)) vertices are 
isolated. 
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Theorem 4.12. (i) If ^ S, then Jl > for every A > 0. 

(ii) If £ S and 1 £ S, then fi > for every A > 0. 

(iii) If £ S and 1^5, then there exists Aq > such that /I = for 
every A < Aq, but fi > for every A > Aq. 

Proof, (i): Trivial, since ^ -E'(A) in this case. 

(ii) : In this case, (psi^) = (p'si^) ~ ^ thus (j4.3p shows that 

-^V5,i(/u;A)^l>0 

as fi ^ 0; hence ipsif^-] ^) is increasing for small and is not a maximum 
point. By Theorem S^^ii) , ^ ^o(A). 

(iii) : By Lemma 14.91 A(/x) tends to oo as — > oo, and to either 1 or oo 
as /i ^ 0. Hence Ai := inf^>o A(/x) > 0. If A < Ai, there is thus no positive 
solution to ()2.5p . so E(X) = {0} and /I = 0. The existence of Aq > Ai as 
asserted now follows from Theorem 14.51 and Theorem 14.61 □ 

In case (iii), /U*(Ao) = by Theorem 14. 6| it is possible both that Eq{Xo) = 
{0} so that /i*(Ao) = 0, and that |i?o(Ao)| > 1 so that Aq € A and /i*(Ao) > 0. 
We can classify these subcases too. 

Theorem 4.13. Suppose that £ S and 1^5. 

(i) If 2 £ S andS ^ S (S = {0,2, s, . . .} with s > 4, or {0, 2} ), then 
/i(A) = for A < Aq = 1 and /i(A) > for A > 1, with fl*{l) = and thus 
1 ^ A and /1(A) \ as A \ 1. 

(ii) If 2 £ S and 3 £ S, or if 2 ^ S (S = {0,2,3, .. .} or {0, s, . . . } 
with s > 3), then there exists Aq > such that /i(A) = for A < Aq and 
/i(A) > for A > Aq, with /i*(Ao) > = J1*{Xq) and thus Aq G A and 
limA\Ao /"(A) > 0. 

Proof, (ii): If 2,3 £ S, then the Taylor series (psifJ-) = 1 + ^/U^ + . . . and 
0^(/i) = ^ + ^/i^ + . . . yield A(/i) = I — . . . for small fi. If 2 ^ S, then 
1/X{fj,) ^ as /i — > by Lemma HTOl In both cases, 1/A is analytic in a 
neighbourhood of and increases on an interval (0, /io), so A decreases there 
and the result follows by Theorems 14.11 ( i) and l4.12T iii). 

(i): Taylor expansions as in the proof of (ii) show that A(/i) = 1 + ^//^ + . . . 
(when 4 G 5) or X{fi) = 1 + ^^'^ + . . . (when 4 ^ 5) so A increases for small 
fj,, say in an interval (0,/io)- By Lemma [4.101 ^5 increases in (0, /io), and 
since V'5 is continuous at we have ipsit^) > "05(0) for < ^ < fiQ. 

However, we also have to consider larger //, and we use Lemma 14.141 
below which implies that Ai := inf^>^j(, A(/i) > 1. It follows that if A < 1, 
then E{X) = {0}, and in particular /i*(l) = 0. Similarly, if 1 < A < Ai, 
then E{X) = {0,/i} for the unique fi £ (0,/io) with A(/x) = A; in this 
case, ipsifJ-) > "05(0) and we have /2 = /x > 0. Finally, by Theorem 14.61 
/2(A) \ /i*(l) = as A \ 1. □ 
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Lemma 4.14. Under the assumptions 0,2 £ S and 1,3 ^ S of Theo- 
rem \4.1^ i), A(/i) > 1 for every fj, > 0. 

Proof. By (j4.6p . the claimed inequality is equivalent to (p'^i^J-) < /^^^(/i), 
where (/'^(/u) = J2k&s /{k — 1)!. First, we use the trivial estimates 

(l)'s{tA < 'Ai>o\{i,3}(/^) = - 1 - 

and 

We may verify numerically (by Maple, or otherwise) that e'^ — 1 — ^/i^ < 
/x(l + ^//^) for < < 3.38, and thus the claim holds in this range. 

For larger n, we split the sum for 4>'g{iJ-) into two parts. With K := [/u^J , 
we have that 

keS,k<K^ ' keS,i<k<K^ ^keS,k>4: 

K 

= fi + —{(t>s{f^)-l-kfJ-^) 

< - ^/^'- (4.7) 

For k > K we use the Chernoff bound for the Poisson distribution, see e.g. 
[l3 . Theorem 2.1 and Remark 2.6]: 

fc— 1 ^ A: 

feScS, fc>X+l ^ ^ k=K 

< exp(/x - log(K/^) + K-fi) = (^^y. (4.8) 

For > 3, K/fi = [fJ.'^l/n > 9/VW > e, so the sum in (j4.8p is less than 1. 
We combine this with (j4.7p to obtain ^^(/i) < ficpsifJ-) for ^ > 3. □ 

Theorem 4.15. T/ie set {A : /I is noi analytic at A} of phase transitions is 
at most countable. Each phase transition is of one of the following types. 

(i) A jump discontinuity: A G A and /i=i,(A) < /I* (A). 

(ii) A continuous phase transition with /I* (A) = /i*(A) = but /I(A') > 
for A' > A. This can happen only at X = 1, where it happens if and 
only if 0,2 gS but 1,3 ^ S. 

(iii) A continuous phase transition with fi > 0; in this case, A = A(/u) for 
some jjL > with A'(^) = but A increasing in a neighbourhood of fi; 
thus fi is an inflection point of X. 

Examples of type (i) and (ii) are given in Section [6l We do not know 
whether (iii) actually occurs. 

Proof. Since A(/I(A)) = A when ju(A) > and A is analytic, the implicit func- 
tion theorem shows that fi is analytic at every point where it is continuous 
and positive and A' (/I) / 0. Hence we have only the three given possibilities; 
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in (iii), A has to be increasing in a neighbourhood of n since otherwise we 

would have a jump discontinuity by Theorem I4.11( i). 

The further characterization in (ii) follows by Theorems 14.121 and 14. 131 
The number of jump discontinuities is countable by Theorem 14.81 and so 

is the number of inflection points of A, while there is at most one phase 

transition of type (ii). □ 

In Case (ii), by the proof of Theorem l4.13l X{fJ.) = l+C;U^+o(/x^) as ;U — > 0, 
so /2(A) ~ c'V A — 1 as A \ 1 and we have a square-root type singularity. 
In Case (iii), provided it happens at all, if the inflection point is then 
A(/i) — A(/Uo) ~ c{iJ, — //o)™" as fi ^ fiQ for some odd m > 3 (presumably 
m = 3), and thus /i(A) — /Uq ~ c'(A — Xq)^/"^ as A ^ Aq = A(/io)- 

Theorems 14.111 and 14.151 show that phase transitions, except the possible 
one of type (ii), occur when A ceases to be increasing at some points, or at 
least almost ceases to be, in the form of an inflection point. (Recall that 
A — > oo, so it increases in the long run.) We have no criterion for when this 
happens, but it seems likely that it occurs whenever there are large gaps in 
S. 

Problems 4.16. Several open problems remain. For example: 

(i) Is there ever any phase transition of type (iii) in Theorem 14.151 (with 
an inflection point of A)? 

(ii) Does the set of phase transitions lack accumulation points? In other 
words, if there is an infinite number of phase transitions (as in Ex- 
ample 16. 9p . can we always order them in an increasing sequence 
Ai < A2 < . . . with A„ — > 00? 

(iii) Is |£'o(-^)| always 1 or 2? In the latter case, is always |-E(A)| = 3, 
as in Example 16.81 with one intermediate point that is a minimum 
rather than a maximum of ips,i and '05,2? 

5. MONOTONICITY 

We begin with a general result that is a simple consequence of standard 
results. Recall that if X and Y are two random variables, we say that X is 
stochastically smaller than Y, and write X <st Y, if F{X > x) < F(Y > x) 
for every real x; it is well-known that this is equivalent to the existence of a 
coupling {X' ,Y') of {X,Y) with X' < Y' a.s. See, for example. Section 
IV. 1]. 

Lemma 5.1. Let Y be a random variable on Z>o with a probability gener- 
ating function (pviz) = J2i>oPi^^ ^^"'^ '^^ finite for all z, and let, for /i > 0, 
Y^ have the conjugate (or tilted) distribution F{Y^ = k) = pkjjf' /(pvilj)- 

(i) /// : Z>o ^ M is a non- decreasing function such that E |/(y^)| < cxd 
for every /U > 0, then ^E/(y^) > 0, with strict inequality except in 
the trivial case when f is constant on {/c : p^. > 0}. 

(ii) // m < 1x2 then Y^^ <st Y^^ . 
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Proof, (i) : 

_ Ek kf{k)pkfi' (Efc /(fe)Pfc/) (Efc fefffc/i^) 

= E(F^/(y^)) -E(y^)E(/(y^)) = Cov(/(y^),y^)) > o, 

since, as is well-known, the two non-decreasing functions /(^) and of 
are positively correlated, for example by a calculation of E [{f{Y^) — 
f{Yp){Y^ — Yp'j > with Y^ an independent copy of y^. The same proof 
yields strict inequality if f{j) ^ f{k) for some j, k with Pj,pk > 0. 

(ii): By (i), P(l^ > j;) is a non-decreasing function of for every x. □ 

Let, for /i > 0, ~ Po(^) and X^^s ~ Pos{fi). Applying Lemma [STT] to 
Xi^s, we obtain the following. 

Theorem 5.2. (i) If f : 7,>q ^ R is a non- decreasing function such that 
E < oo for every /U > 0, then 

^^f{x,,s) = ^Wi^,) I x,es)> 0, 

with strict inequality except in the trivial case when f is constant on S. 
(ii) If m < fi2 then X^^^^s <st ^/Z2,5- 

This shows, in conjunction with Theorem l2.1l that the asymptotic degree 
distribution of Gn^\„/n]S is stochastically increasing in /1(A) and thus, by 
Theorem 14.51 in A. In particular, the asymptotic edge density, which is 
given by EXp;(^x),s = ^if'-W)^ is increasing function of A (except that it 
is constant in the trivial case \S\ = 1). This holds for finite n too. 

Theorem 5.3. If < Pi <P2 < ^, then e{Gn,pi;s) <st e(G„,p2;-s)- 

Proof. This is another application of Lemma 15.11 since it follows from (jl.ip 
that e{Gn,p-s) has the conjugate distribution with Y = e(G^ i.^). 

□ 

Unfortunately, if we consider the entire random graph Gn^p-s (and not 
just the number of its edges), it is in general not stochastically increasing 
in p. 

Example 5.4. Let n = 4 and let S = {0, 2} (or the set of all even numbers). 
The 5-graphs are, ignoring the labelling: (i) S4, the empty graph with no 
edges, (ii) Gs + Ei, a 3-cycle plus an isolated vertex, (iii) C4, a 4-cycle. 
We have F{Gn,p-s = E^) ^ 1 as p ^ and W{Gn,p-s = C4) ^ 1 as p ^ 1. 
Hence, if /(G) is the number of 3-cycles in G, then E f{Gn,p;s) = ^{Gn,p;S = 
Cs + El) tends to both as p — > and p ^ 1, so this expectation is not 
monotone in p. 
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Problem 5.5. Is the random multigraph G* ^.g defined in Section [7| stochas- 
tically increasing in ul (Its number of edges is, by the same argument as 
for Theorem 15.31 ) 

For the existence of a giant component, we note that the crucial quantity 
Q(/u) in (13. 2p is not always monotone in not even in the classical case 
S = Z>o (when Q{p.) = — /i). Nevertheless, the condition Q{^) > is 
monotone. 

Theorem 5.6. If fii < ^2 and Q{^i) > 0, then Q(/U2) > 0. 

Moreover, assuming S ^ {0,2}, if Xi < A2 and thus /i(Ai) < /u(A2), 
i/jen ^(/2(Ai)) > ^(/I(A2)) and 7(/i(Ai)) < j{Jl{X2)). Hence, if Gn,Xi/n;S 
has a giant component, then so has Gn^\2/n;S f^'^ -^2 > Ai, and it is 
(asymptotically) at least as large. 

Proof. The condition Q{^) > is equivalent to /i0^(/i)/0^(/i) > 1. Since 
05 (^) = (t^S^lilA^ 



l^fPsil^) _ ^'^5-l(/^) 



IEX^,5-i, (5.1) 



which is non-decreasing by Theorem 15. 2 f i). 

The monotonicity of ^ and 7 follows from the branching process inter- 
pretations in Remark 13.51 together with the stochastic monotonicity Theo- 
rem [521^11) (for both S and 5 — 1) and Theorem 14.51 □ 

Remark 5.7. Theorem 13.11 and (j5.ip yield the curious relation that, pro- 
vided S 2 {0,2}, Gn^x/n-s a giant component if and only if i/5_i(/I) := 
KX^^g^i > 1, with fi = /i(A) calculated for S. Cf. Remarks 13.51 and 13.71 

It follows similarly from Remark [321 that the existence and size of a /c-core, 
for any fixed /c > 3, is monotone in A. 

6. Examples 

Example 6.1. S = Z>o, the Erdos-Renyi random graph. We have in 
this much studied case that (psif^) = The characteristic equation ()2.5p 
becomes fi = ff^/X, with solutions fi = 0, ^ = X. By (12. 7p . ipsifJ') = ^/^i so 
V'5(0) < i^siX). Therefore, and in accordance with Theorem 14. 12( 11) . /x = A, 
so that the number ni of vertices of degree i satisfies Ui/n A*e~'^/i!, i > 0. 
This is a simple instance of Theorem l4. 1 1 ( ii) . There is no phase transition of 
Gn,K/n;S- We have ^5,i(/^; ^) = /"-/"^/(2A) and ^psAl^] A) = ^^(log(A//i) + 
1) . Details of th e ap plication of Theorem 13.11 to this well understood case 
may be found in [l5|]. Similarly, the application of Theorem 13.61 is described 
in 0. 

Example 6.2. S = {s}, where s > 1. We have (j^sin) = /i*/s! and the 
characteristic equation ()2.5p becomes s = /i^/A, with solution fi = y/sX. 
However, in this case, the value of /I is in fact immaterial, since Po5(/i) 
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is a point mass at s for every > 0. Moreover, the graph Gn^x„/n;S is a 
random regular graph with all vertices of degree s. It is immediate that 
= s(s — 2), and so there exists a giant component if s > 2, and not if 
s = 1. In fact, if s = 1, the graph consists of isolated edges only, while if 
s > 3, it is well-known that the graph is connected with probability tending 
to 1, see Bollobas 0, Section VII.6]. In the remaining case s = 2, the 
graph consists of cycles, of which the largest has a length that divided by 
n converges to some non-degenerate distribution on [0, 1], see, e.g., Arratia, 
Barbour and Tavare [H; this is thus an exceptional case where we do not 
have convergence in probability of the proportion of vertices in the giant 
cluster, as in Theorem 13.11 

Example 6.3. S = 2Z>o, the even numbers. In this case, 

'^'S(/^) = J^^ = cosh/i. 

The characteristic equation (|2.5|) is 

sinh fi fj? 
cosh fi A ' 

so either /i = or 

A = A(/z) = -i^. (6.1) 
tanh// 

Since A(;u) = /u/ tanh/u increases (strictly) from 1 to oo for /i E [0,oo), it 
follows that: if A < 1, ;U = is the only solution, while if A > 1, there is also 
a positive solution. We have 

"05 (a*) = log (cosh /i) — i^tanh//. 

Therefore, 

sinh(2/x) — 2/i 
4 cosh^ fi 

for /i > 0. Hence, "05 (a*) > V'5(0) ^or /i > 0, whence fl is the unique 
positive solution of (16. ip when A > 1, cf. Lemma HTTUl Theorem 14. 1 l( ii) and 
Theorem I4.13( i) . 

We thus have a continuous phase transition at A = 1 with /2(1) = 0; there 
is a unique fl (and thus Theorem 12.11 applies) for every A, and /2(A) is a 
continuous function, but it is not differentiable at A = 1. This is the only 
phase transition, and A = 0. 

The asymptotic edge density (i.e., the number of edges per vertex, see 
dZH])) is 

u{'fl) = /itanh/I. (6.2) 

Since 1^5, Theorem 13.11 shows that there is a giant component as soon 
as /2 > 0, i.e., if A > 1. In fact, it is easily seen that 

\tanh n 



, , , , smni zi/i I — zu, 



(5(/i) = /itanh|U 
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One may study a random even subgraph of a general graph G. It turns out 
that the random even subgraph with parameter p G [0, |] is related to the 
random-cluster model on G with edge-parameter 2p and cluster-weighting 
factor q = 2. When G is a planar graph, the random even subgraph may 
be identified as the dual graph of the +/— boundary of the Ising model on 
the (Whitney) dual graph of G with an appropriate parameter- value. This 
relationship is especially fruitful when G is part of a planar lattice such as 
the square lattice 7?. See [fi] for a general account of the random-cluster 
model, and 0] for its relationship with the random even subgraph and the 
Ising model. 

Example 6.4. S = {1,3,5,...}, the odd numbers. This time, (j^sif^) = 
sinh /i, and the characteristic equation is 

// cosh n fj? 
sinh fj, A 

with A(/i) = /itanhyu. Since 0^5 and A is increasing, the unique solution 
/I is given as the unique positive solution of /utanh/i = A. Cf. Theorems 
I4.12r i) and l4.1lT ii). There is no phase transition. This time, 

^(^) = , (/itanh^- 1). 
tanh fi 

Thus Q{fi) > if and only if /I tanh /2 > 1; since /z tanh /I = A, it follows 
that there is a giant component for A > 1, and not for A < 1. (This also 
follows from Remark 15.71 and (j6.2p .) In the critical case A = 1 we have 
/itanh/z = 1 and numerically fl w 1.19968 and (asymptotic) edge density 
„{f2) = ft^/X = J1^ ^ 1.43923. 

Example 6.5. S = {1,2,3,...} = Z>i, graphs without isolated vertices. 
We have that 4>s{^i) = — 1, and the characteristic equation is 



- I A 

Since ^ 5, we seek strictly positive solutions, which is to say that A = 
A(/i) = fi{l — e~^). Since fi{l — e~^) is increasing on (0, oo), there is a unique 
such solution /I for every A > 0. Cf. Theorem 14.11 ( ii) . 
We have that 

— 1 

Thus Q{fi) > if and only if /I > 1, which is to say that A > A(l) = 1 — e~^. 
There is a giant component when A > 1 — e~^, and not when A < 1 — e~^. 
In the critical case A = 1 — e~^, Ji = 1 and the critical (asymptotic) edge 
density is i^(/2) = /IVA = e/(e - 1) « 1.58198. 

Example 6.6. S = {0,1}, matchings. We have that (psifJ-) = 1 + A*) and 
the characteristic equation is 

1 + /X A ■ 
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Either = or A = A(/i) = ^(1 + /i), so the solutions for given A are /u = 
and ;U = — i + -i/A + i. Since A is increasing, Theorem I4.11( ii) apphes and 



shows that fl = ~^ + y ^ ~^ \ for all A > 0. This can also easily be verified 
directly, using 



2(1 + /.) 

which yields ip'si^J^) > for fi > (cf. Lemma l4.10p . 
By Theorem 12.11 as n ^ oo, 




n (t)s{V) 1 + /I 

Obviously there is no giant component. Indeed, Q(/x) = — Po5(/x){l} < 
for /i > 0. 

Example 6.7. S = {0,2}, isolated cycles. We have that (psilA = 1 + ^/U^, 
and the characteristic equation is 



1 + A • 

Therefore, either = or A = A(/i) = 1 + ^/i^, so that the solutions for 
a given A are = and, when A > 1, /x = ^J2{\ — 1). Again, A is an 
increasing function, and so is 

Ml^) = log(l + 2^ 



l + i/x2' 

by Lemma 14.101 or direct calculations. Thus, see Theorem 14.11 f ii) . 

when A < 1, 



Y^2(A - 1) when A > 1. 



We thus have a continuous phase transition at A = 1, of the same type as in 
Example 16. 3^ see Theorem 14. 12r iii) and Theorem I4.13l fi). There is no other 
phase transition. 

It is easily seen that = for all /x, which may be interpreted as 

saying that the random graph is, in a certain sense, critical whenever A > 1. 
If we remove the isolated vertices, and condition on the number of remaining 
vertices, we obtain a random regular graph with degree 2. Hence, for /2 > 0, 
i.e., for A > 1, we see that the largest component behaves as for S = {2}, 
see Example 16.21 with convergence of v(r„ to a distribution but 

not to a constant. 

Example 6.8. S = {0, 3}. This time, (f^sifJ^) = 1 + g/f^^i and the character- 
istic equation is 

l + i//3 A- 
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Either ^ = 0, or 

A = A(^) = (6.3) 

2^ 

This is a convex function of ^ with a minimum of 3^/^ at the point /i = 3^/^. 
Hence, the characteristic equation has no positive root when A < 3^/^, one 
such root if A = 3^/^, and two such roots if A > 3^/^. We have 

i + g/X 

Unhke the previous examples, "05 is not monotone, cf. Lemma [4. 101 In fact, 

i^si^^h = log(l + i) - i < = V5(0), 

so the correct root is /I = (rather than 3^'''^) when A = 3^^^. The function 
"i/i has a minimum at ^ = 3^/^; V decreases on [0,3^/^] and increases on 
[3^/^,oo). There exists thus a unique /io > 3-*^/^ such that i^sifJ-o) = 0) and 
we set Ao = 2(1 + i/iQ)/;Uo. (Numerically, ^uq ~ 2.03134 and Aq = A(/io) ~ 
2.36002.) We deduce that /i = for A < Aq while, for A > Aq, /I is the 
largest root of (16. 3p . For A = Aq, there are two roots /x of ()6.3p with the 
same value of -05(^)1 so we have a jump phase transition and Theorem 12.11 
does not apply; see Theorems I4.12( iii) and l4.13T ii). There is no other phase 
transition, and A = {Aq}. 

Since 1 ^ 5, by Theorem 13. l( ii). there exists a giant component whenever 
A > Aq. Indeed, 

l + i^3 

for every > 0. 

Example 6.9. 5 = {1, 2, 4, 8, . . . } = {2^ : j > 0}. We claim that as j oo, 

A(2^x) = 2^x2(1 + o(l)), (6.4) 

for every x £ (2/e,4/e). (In fact, this holds uniformly on every closed 
subinterval of (2/e, 4/e).) It follows that if a = 4/e and e > is small 
and fixed, then for large j, A((a — e)2^) « (a — e)^2-' and A((a + e)2^) = 
A((a/2 + e/2)2J+i) (a + e)^2^-'^, so A drops by a factor of about 2 in 
the vicinity of a2^ . Consequently, for all large j, there is an interval Ij C 
{^{a — e)2^ , (a + e)2-') where A is decreasing, and thus by Theorem 14.111 there 
exists A G A such that /i*(A) > max/j > 2^ . Hence the set {/2*(A) : A G A} is 
unbounded and thus infinite, so A is infinite and there is an infinite number 
of phase transitions. 

To verify ()6.4p we show that if = 2^x G (02-^"^, 02-^) then (j)s{fJ-) = 
X^fce^A*^/^' and ^^(/i) = J2kes ^l-'-^ / ~ dominated by the terms 

with k = 2^: 

M^^) = (l + o(l))^, <A5(/^) = + (6-5) 
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as j — > oo; this yields that ~ {2^ / fi)(j)s{n) and 

Finally, to show (j6.5p . we observe by Stirling's inequality that, as A; ^ cxd, 
and thus, with ki = 2'', 

0/|;:Hi-a,).-(|^)-(f)-" 

which for ^ = 2^x £ {a2^~^ , a2^) is exponentially small if i > j and expo- 
nentially large if i < j. The estimate ()6.5p for (psil^) follows, and a similar 
calculation with /cj = 2* — 1 yields the result for (j)'g(fi). 



7. MULTIGRAPHS 

As explained at the end of Section [21 we shall count multigraphs with 
certain properties, and shall later relate our conclusions to simple graphs. 
Let ^* be the (infinite) set of all multigraphs on the vertex set {1,2, ... ,n}, 
and let be the subset of 5-multigraphs on {1,2, ...,n} (we extend 
the definitions above to multigraphs in the obvious way, noting that a loop 
counts two towards the degree of the vertex in question) . 

Let > 0. We define a random multigraph G* ,^ by taking Po(z^) edges 
between each pair of vertices and Po(^i^) loops at each vertex, these random 
numbers being independent of one another. It is easily seen that this is 
equivalent to assigning to each multigraph G € ^* the probability 

P(G;_^ = G)= w{G)u<^'>e-'''''/\ (7.1) 

where 

w{G) := 2-'llj\-^^, 

i>2 

with i the number of loops of G, and rrij the number of j-fold multiple edges 
(including multiple loops). That is, rrij = aj + bj where aj is the number 
of distinct pairs of vertices joined by exactly j parallel edges, and bj is the 
number of vertices having exactly j loops. See, e.g., Janson, Knuth, Luczak 
and Pittel 0. 

Note that the total number of edges e(G* j^) is Poisson-distributed with 
parameter {^v + ^nv = ^n'^i'. We further define the random 5-multigraph 
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G* j^.^ by conditioning G* ^ on being an 5-multigraph. Thus, for any multi- 
graph G G gis, by (mD, ' 

= G)= -^w{G)v<^\ (7.2) 

ri,u;S 

where 

We shall assume, of course, that ^*.^ 7^ 0. It is easy to see that this holds 
for all n if 5 contains some even number, but if all elements of S are odd, 
then n has to be even. We tacitly assume this in the sequel. 

If the multigraph G G ^* is simple, i.e. has no loops and no multiple edges, 
then w{G) = 1 and yields P(G;^^ = G) x x P(G„,p = G) when 

f = p/(l — p), i.e., p = u/{l + v). Hence, assuming this relation between v 
and p, G* conditioned on being simple has the same distribution as Gn,p- 
Conditioning further on being 5-graphs, we obtain the following. 

Lemma 7.1. // v = p/{l —p), then 

Gn,p;S = {G*n,u;S I ^l^^.g simple) . 

We are interested in the case np ^ X < 00, and note that np ^ X and 
nv ^ X are equivalent. 

We shall also use the configuration model for random multigraphs with 
given vertex degrees introduced by Bollobas [4| , see Bollobas d. Section II.4]. 
(See Bender and Canfield [3] and Wormald [Ig, [l9] for related arguments.) 
To be precise, let us fix the vertex degrees to be some non-negative integers 
di,d2, ■ ■ ■ ,dn (assuming tacitly that is even); equivalently, we fix a 

degree sequence d = (di)". We attach di half- edges (or stubs) to vertex i. 
The total number of half-edges is thus 2 A'" := Y17=i '^i-> ^"^^ ^ configuration is 
one of the (2iV-l)!! = (2iV)!/(2^iV!) partitions of the set of half-edges into 
N pairs. Each configuration defines a multigraph in ^* by combining each 
pair of half-edges to an edge; this multigraph has vertex degrees di, . . . ,dn 
and N = ^ di edges. By taking a uniformly random configuration we 
thus obtain a random multigraph G* ^ with the given degree sequence d. 

It is easily seen that every multigraph G € with the given vertex de- 
grees di, . . . ,dn arises from exactly w{G) YVi=i ^i' configurations. We obtain 
therefore that the contribution to Z* ^^.g in (17. 3p from a set of multigraphs 
with given vertex degrees di, . . . , d„ G 5 is given by summing / HiLi 
over all corresponding configurations. In particular, since the number of 
configurations is (2A^ — 1)!!, the contribution to Z* ^.^ from all multigraphs 
with vertex degrees di, . . . ,dn G S equals 

(^^-1)"-^ (7 4) 
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Moreover, for given d = ((ij)", the factor / YYi=i ^i^- ^ constant, so 
by (17.ip , the probabihty that G* belongs to any given set of multigraphs 
with this degree sequence d is proportional to the number of corresponding 
configurations. Consequently, if di{G) denotes the degree of vertex z in a 
(multi)graph G, and d(G) := {di{G))f^i for G G we obtain the following 
well-known fact. (This is another reason for the weights w{G) in (jT.ip .) 

Lemma 7.2. For any given degree sequence d = (dj)" and any v > 0, the 
random multigraph G* conditioned on having degree sequence d has the 
distribution given by the configuration model; in other words, 

{Gn,u I ^iGn,u) = d) = G* 
As a consequence, if every di G S, the same holds for G* ^.g. 

We are interested in the case v = A„/n, with A„ ^ A > 0. 

Theorem 7.3. The results of Theorem \2. l^ i)^( Hi) hold with Gn^\^/n;S ^6- 
placed by C^An/n cS- Furthermore, 

i logP(G;,„/„ is an 5-graph) ^ ^^(/S) - (7.5) 
and n-MogZ*^^^/„.5 ^ i;s{J2). 

We will prove Theorem 17.31 in the following section, and then obtain The- 
orem 12.11 as a consequence using Lemma 17.11 and the following technical 
result. 

Lemma 7.4. // A„ — > A > 0, then liminf„_+oo IP(G* ^ ^^.^ is simple) > 0. 

8. Proof of Theorem 17.31 

For notational convenience, we shall consider only the case A^ = A for 
all n, while noting that our estimates may be extended to the general case 
Xn — > A. (The "constants" below generally depend on A, but they may be 
chosen uniformly for A lying in any compact subset of (0,oo). Uniformity 
as A —> is less obvious, and perhaps not always true, but it is remarked a 
few times when it is important for later proofs.) 

Let A/^ denote the set of all n = (no, ni, . . . ) G Ms such that: Ylj = ^ 
and Ylj is even. We write z{n) for the contribution to ^^^.^ from all 
multigraphs with nj vertices of degree j, which is to say that 

K,X/n;S = E ^(") (8.1) 

and 

n^iG:,X/n;s) = n) = n e Ml (8.2) 

n,X/n]S 
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By ([73]) with N = ^Ylj J^j and v = X/n, 



n^- {Ejjnj - 1)!! ^A^^ 



We note by Stirling's formula that 

and it is easily verified that 

(2iV -!)!!> i^—j , N>0. (8.5) 

Let n = n(n) be a mode of the random sequence n(G* x/n-s^^ (!8-2p . 
a sequence in A/"^ that maximizes z(n). (In the case of a tied maximum we 
make an arbitrary choice.) We write N = N{n) := ^ '^jj^j- 

We begin with a coarse but useful quantitative estimate, obtained by 
considering only regular multigraphs. Let := (5js)i^0' '^here 5is is the 
Kronecker delta. In the following lemma we take an even number s G 5, if 
S contains such a number. If not, we pick an odd s G 5 and must then, as 
noted in the introduction, restrict ourselves to even values of n. 

Lemma 8.1. Let s G S, and assume that s is even if possible. Then 

C5(s)" < z{nes) < z{n) < Z;,/^^^ < (8.6) 

As a consequence, 

^iG*n,x/n is an 5-multigraph) > c^, 
and, for any set 7i of multigraphs, 

Proof By (gS]) and (fH^ll . 

■\\ns/2 f rm \\'n-s/2 

-) >(s!)-"( — .-^ 

which yields the first inequality in (|8.6|) with C5 = (s!)~^(sA/e)*/^. The 
second and third inequalities are trivial, and the fourth follows from (|7.3p 
(with u = X/n), which yields Z* ^^^.^ < e"'^/^ ^ 
By ^ and ([SSD, 

IP(G:,A/n is an 5-multigraph) = e-'"/^ Z^x/n^s > 

with cg = cse"'^/^. Consequently, (j8.7p follows with C7 = Cg ""^ by the defini- 
tion of conditional probabilities. □ 

Lemma 8.2. There exists a constant B = B{S, X) such that 

z{n) < e-z(a) < e-"Z;,/„^^, 



z{nes) = (s!)-"(ns - 1)!! ( -Y'^^ > (s!)-"f— • - 

\n/ \ e n 
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where N = ^ Yli^ni. Hence, IP(e(G; ^^z^.^) > Bn) < e"". 
More generally, for any x > Bn, 

Moreover, for any Aq > 0, the constant B can be chosen uniformly for all 
A< Ao. 

Proof By (I7T1) . 

5; z{n)< ^(n)=e"'^'lP(KG:,,/J>x). (8.8) 

n€Afg:N>x neAf":N>x 

Since the number of edges e(G* ^^^) ~ Po(^An), it follows by standard 
Chernoff estimates for the Poisson distribution, see e.g. jl2! . Corollary 2.4 
and Remark 2.6], that if S > 4A and x > Bn > 4An, then 

IP(e(G';,/J >x)= P(Po(iAn) > x) < e'^ (8.9) 

We choose 

B > max{4A, ^A + 1 - log 05(5)}, 
and find by ([M]), <n and (US]), since {B - l)n < {B - l)x/B = x- x/B, 
z(n) < e"V2-x-nlogc5(s)^^g^ < e(^^i)""^z(S) < e^^'/^z(fi). 

The results follow by this and (j8.2p . □ 

Let n G Ms, and let j and /c be two different indices in 5 such that > 2, 
and define n' G A/^ by n'- = Uj + 2, n'^ = — 2, and n'j = nj for i 7^ j, k; in 
other words, we replace two vertices of degree k by vertices of degree j. By 
(I83D . with = ^ Ei and A^' = ^ Ei = iV + j - fc, 

z(n') _ (2iV' - 1)!! Ujlnkl k9 (\\N'-n 



z{n) (2iV-l)!! nMn'^l j^? Vn. 

For n = n and any j,k ^ S with > 2, this quotient is < 1. Hence, for all 
j,k ^ S (also, trivially, if < 2 or j = k), 

Mnk-m' < (% + l)(% + 2)i!2 (^y~' (1 + 0(|j - A:|ViV)) . (8.10) 

Furthermore, in the case k > j, we have N' < N and 

{2N - 1)!! < {2N' - 1)!! (2iV)'=-^ 
and we obtain in the same way the sharper inequality 

Mnt - < (%■ + !)(% + 2)j!2 (^)'~'. (8.11) 
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By Lemma l8.2[ 

N < Bn. (8.12) 

Now let n — > oo. Since N/n is bounded by (j8.12p . each subsequence has a 
subsequence such that N/n converges. Consider such a subsequence, and 
assume that 2N /n v >0. Furthermore, let pj := nj/n. Then p := {pj)'^' 
is a probability distribution (the distribution of the degree of a random 
vertex in a graph G with n{G) = n). Since the mean of this distribution 
is YljjPj — 2A^/n, which is bounded by (|8.12p as n — » oo, this sequence of 
distributions is tight, and by taking a further subsequence we may assume 
that the distributions converge, i.e. that rij/n — > pj for some probability 
distribution (pj)^ and every j > 0. Clearly, this probability distribution is 
supported on S in that pj = when j ^ S. 

We treat the cases u > and 1^ = separately. Assume first that z/ > 0. 
Divide (KIM by and let n ^ oo to find that 

plk\''<p]j\HX,.)'^-^, j,kGS, 

and thus 

plk\\Xu)->^ < p^jl^ {Xur^, j, k G S. (8.13) 

Interchanging j and k we obtain equality in (j8.13p . Writing C| for the 
common value, and ^ := ^fXu, we deduce that 

Pj = Csj. j^S. (8.14) 

If, instead, v = then X]j>i^j — — o('^)) so no/n — > 1 and nj/n 0, 
i > 1; hence pg = 1, and Pj = for j > 0. (Thus, = implies that G S.) 
Equation (j8.14p holds in this case also, this time with /i = 0. 

Hence, (j8.14p holds in all cases. Summing over j and recalling (j2.2p . we 
find that 1 = Csfpsil^) and thus Cg = (psilJ-)"^- In particular, (psifJ-) > 0, 
another demonstration that /i = is possible only when £ S. 

In summary, along the selected subsequence, 

p, := ^ - p, = ^ = Po5(/x){j}, J G 5, (8.15) 

which is to say that every subsequence possesses a subsequence along which 

p = ip,)^ ^ Po5(m) (8.16) 

for some fi. We next identify /x, and show that it is the same for all subse- 
quences. 

We constructed n' above by changing by 2 the degrees of two vertices; 
the reason was that this ensures that ^ • in'- remains even. If j and k have 
the same parity, i.e. j — k = (mod 2), then we may also argue as above 
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changing just one vertex degree from fc to j. If further k > j, this leads as 
in (jS.lip to the inequaUty 

ukkl < (% + l)j! . (8.17) 

We consider again a subsequence along which (j8.16p holds for some /x. For 
every k we apply (j8.17p with j the smallest number in S of the same par- 
ity as k. Using (|8.15p for these (at most two) j, we obtain, with /i„ := 
(2AiV/n)V2 ^ (Azy)V2 ^ uniformly for all k£S, 

n. < l^n{p^+ oil)) jlf^t' = ^;^Mn-^(/^^' + 
since < /i + 1 for large n. Consequently, for every exponent r > 0, 

^ Fp, = E ^^TT - + o(i))05(^)-i E + 1)' = 0(1)- 

fce5 fee5 fce5 

In other words, for every r G (0,oo), the distributions p have rth moments 
that are uniformly bounded in n. It follows that all moments converge in 
()8.16p . i.e., for every r > 0, 

^ Fpfc = E ^'■^r ^ E ''^Pk = E Po5(m){^}- (8-18) 

kes kes kes k&s 



In particular, r = 1 yields, using ([27 

— = E^^-E^P. = ^#r- (8-19) 

On the other hand, we have assumed 2N/n v and /i = -v/Aj^, whence 
I' = /i^/ A, so we have the consistency relation 

In other words, with E{X) defined by (|2.6p . we have [i £ E{X). 

We summarize the result so far. Each subsequence of n possesses a sub- 
sequence such that (|8.16p and ()8.18p hold for some ^ G E{\), i.e., /i > 
satisfies (I8.20p and, further, fi = only if £ S. 

The next step is to find the right solution of ()8.20p in the case when E{X) 
contains two or more points. 

We continue to consider a subsequence for which (|8.16p holds. By apply- 
ing again ()8.17p with j the smallest odd or even number in S as appropriate, 
and using rij < n and ()8.12p . we see that for some constants Cg, Ciq, 

nkk\ < CgnC^Q. (8.21) 
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If /c > logn, then by Stirling's formula, for large n, 

log(A;! C"'^) > A; log A; - k{l + log do) > 2k > log{Cgn), 
and thus (j8.2ip yields < 1. Consequently, for large n, 

nfc = for all k > log n. (8.22) 
Let us now estimate z{n) = maxn2;(n). By (|8.3p and Stirling's formula. 



recalling ([87 

logz(n) = nlogn — n + 0(log?7,) — ^^(f^i logfij — nj + 0(log(nj + 1))) 

i 

-^Ui \og{i\) + iV(log(27V) - 1) + 0(1) + N log(A/n). 

i 

By ()8.22p . we only have to sum over i < logn, and thus the sum of all O 
terms is O(log^n). Thus, with y := "^^ipi = 2N/n, 

— logz(n) = logn — 1 — y^pjflogn + logpi — 1 + logi!) 
n ^-^ 

+ |(logy + logA-l) +o(l) 
= -J^paog(p,i!) + |(logy + logA-l) +o(l). (8.23) 

i 

For each i, pi ^ Pi by (l8J6]) . and in addition, by ([8J9l) and ([H^O]), 

2/ = J^^K^J^^ft = /"VA. (8.24) 

i i 

Now, xlogx > — e""*^ on [0, oo), whence log (pji!) > —e~^/i\, and by (|8.2ip . 

Pi\og{pii\) < Pi log i! < piilogi < pii^ = 0(i^Cio/«!) • 
Consequently, by dominated convergence, 

- log z(fi) ^ - J]lJilog(p,i!) + ^ (log ^ + log A - l) . (8.25) 

i 

Furthermore, by (f8?T5]) and (fOil) . if /i > 0, 



^Pjlog(pji!) = ^Pj(ilog^ - log05(/i)) = y log/i - log05(/i); (8.26) 
if = this holds trivially with all terms zero. Hence, ()8.25p yields, using 

1 u? 

-log2;(n) — -log/x + log05(/i) + — (21og/i- l) 

n A zA 

= log *.(M)-g = log (8.27) 
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Conversely, take any finite sequence {pi)fLQ with > 0, X^jK = 1 ^'^^ 
Pi = when i ^ S. Define rii = nxi, rounded up or down to integers, 
preserving ^ • rii = n and possibly adjusting two of them by ±1 so that 
irii is even. As n ^ oo, we then obtain as in (18.23p - (j8.25p (but simpler, 
since the sums are finite), with u := Y^- ipi, 

- log 2;(n) ^ - Pilog{pii\) + ^ flog + log A - 1 

i 

Since, by definition, z(ji) is maximal, z{n) > z{n), and thus 

liminf i^^ifM > _ pilog{piil) + ^ flog + log A - iV (8.28) 

i 

We have shown (|8.28p for any probability distribution (p,) on S with finite 
support. More generally, let (pi) be a probability distribution supported 
on S with u := ^^ipi < oo and Pilog{pii\) < oo. For M > min5, 

let pf^) := Pi/Zj<MPj for i < M and apply ^E!) to (pf^^)f£o- I* is 
easily seen that the right hand side of ()8.28p converges as M — > oo to the 
corresponding value for (pi), showing that ()8.28p holds for (pi) also. 

In particular, for any /j, E E{X), we can use ()8.28p with pi = Po5(/i){i} 
given by ([231) and, by ([22 



2 



(t>s{^J) A 

Hence, (l8:28D yields, by the calculations in (lOHD and ([8211), 

liminf > log0,(,) - = V^,(;.), (8.29) 

for every /i G E{X). Comparing this to ()8.27p . we see that if ()8.16p holds for 
some subsequence and some /i G E{X), then this fj, must maximize ipsifJ-) = 
V'5,i(/^jA) over E{X), in other words, /U = /I as defined in Theorem 12.11 In 
particular, this shows that every subsequence possesses a subsequence such 
that (j8.16p holds with a fixed = Jl; hence (I8.16P holds for the full sequence 
of n — > oo, and 

p = {Pj)^ ^Posifl). (8.30) 

Remark 8.3. We have for simplicity considered /i G E{X) only in (I8.29P : 
for general /U > the lower bound obtained from this argument takes the 
form, with v = fi(l)'g{fi)/(j)s{fi), 

log (j)s{n) + ^ (log ^ - = 7ps,2{l^; A) 

defined and studied in Section HI 

We have so far studied the mode p of the degree distribution. We now 
show that the distribution is concentrated close to the mode. 
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Lemma 8.4. For every e > 0, there exists cy = C'j{e) > such that, if n 
is large enough then for every n E J\f^ with dxv (n/n, Po5(/i)) > e, z{n) < 

We will first show a weaker statement. 

Lemma 8.5. For every e > 0, there exists = cs{£) > such that, if n 
is large enough then for every n G Af^ with cIty {'n/n,Fos{Jl)) > e, either 
z{n) < e~'^*"z(n), or there exists n' G with (iTv(ii/ra, n'/n) < 2/n and 
zin) < (1-C8)z(n'). 

Proof. Suppose this fails. Then there exists e > and a sequence n = n^") G 
with n ^ oo, such that (Ity [n/n,Fos(jl)) > e, z{n) = e~°^"^2;(n) and 
z(n') < (1 +o(l))2;(n) for all n' E with dTv(n/n, n'/ri) < 2/n, i.e., for 
ah n' G Mg with jn, - n-| < 4. 

We now repeat much of the arguments presented above for the mode n. 
First, we obtain that (fSTOj) . KU\} and (fSTf]) hold for these n, with an 
extra factor 1 + o(l) on the right hand sides, uniformly in all j,k ^ S (with 
k > j for (jS.lip and k > j and k = j (mod 2) for ()8.17p ). Furthermore, 
by Lemma 18.21 and the assumption -z(n) = e~°*-"^2;(n), := ^ Si^'^i ^ Bn 
(for large n). It follows, as above, that by considering a subsequence we 
may assume that 2N/n u € [0,oo) and n/n — > p for some probability 
distribution p, where, again as above, necessarily p = Fosil-i') ^r some 
/i G ^(A) and (KTl\\ holds for n. 

Since logz(n) = o(n) + log 2;(n), this shows that ipsifJ') = "^siJj-), and thus 
// = /i, since Jl is assumed to be a unique maximum point. Consequently, 
n/n — > Fosifi), which contradicts cZtv (n/n, Po5(/i)) > e. □ 

Proof of Lemma \8.4\ By Lemma [8.51 if is large enough, dTv{'n-/n, Fos{'p)^ > 
e, and z{n) > e~^^^z{n), there exists n^^^ = n' such that (iTv(ii/n, n^^^/n) < 
2/n and z{n) < (1 — C8)z{n^^^); in particular, z{n^^'>) > z{n). 

If also (iTv(n''^V^) PO'S(/^)) ^ we iterate and find n^^^, and so on. 
This gives a sequence n^*^) = n, n*^"'^^ , • • • , n^^^ , where for I < L we have 
(iTv(n('V^!n'-^"^^V") < 2/n and z(n(')) < (1 - C8)z(n('+i)), while 

dTv(n(^Vn,Po5(/i)) <e. 

If further (iTv(n/n, Po5(/2)) > 2e, it follows that dTv(n/n, n*^^)/n) > e, 
and thus the number of steps L > en/2. Consequently, 

z{n) < (1 - C8)^z(n(^)) < (1 - cs)^z{n) < exp(-ic8en)z(a). 

This proves Lemma 18.41 for 2e, with C7(2e) = min{l, ^e}cs{s). □ 

We now complete the proof of Theorem 17. 3[ Let e > and let, with B as 
in Lemma 18.21 and N = ^ Yli ^^j, 

Ai := {n e MS ■■ N > Bn}, 

Al := {n e Ms : N < Bn and dTy{n/n,Fos{Jl)) > e}. 
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Lemma 8.6. For every e > 0, < = e°("). 

Proof. Suppose n = (nj)j G For each i, < n, and thus the number of 

choices of {ni)^^^ is at most = exp(0(\/nlogn)). Furthermore, 

iui = 2N < 2Bn, and thus = for i > 2Bn and 

SO (^^^)Vn<^<2_BAf be described by a sequence of at most 2By/n numbers 
in the range [y/n,2Bn] (the degrees of the corresponding vertices). Hence, 
the number of choices of {ni)^^^ is at most {2Brif'^^ = exp(0(-y/nlogn)). 
Combining the two parts, |^2l = exp(0(-y/nlogn)). □ 

Now, fix e > 0. By Lemmas E21 E31 and E3 
P (dTv(vr(G„,,„/„^5),Po5(/2)) > e) < e"" + ^ < e~^^\ 

n,\/n:S 

for some ci > and all large n. This proves (12.12P and hence (|2.9p (for 

^n,A/n;5)- 

A similar calculation with e = yields 

and thus logZ* ^^^.^ = logz(n) + o(n), which together with (18.27P implies 
n-MogZ* ;^/^.^ ipsCfi)- By dLSD, this further yields fTHI) . 

Lemma 8.7. Uniformly in all k >0, 

Enfc(G;,/„^5)<Cnne-^«^ 

Moreover, for any Aq, this holds uniformly in X < Xq. 

Proof. Let j be the smallest element of 5 with the same parity as k. Given 
any n G M^, let n' G be given by n'j := + 1, n'^ := — 1 and n',- := nj, 
i j,k (assuming j < k and > 1; otherwise n' = n). By (j8.3p . cf. the 
argument yielding ()8.17p . 

and thus 

Ci2n/2AiVx(fc-j)/2 

nfcz(n) < — — 2;(n). (8.31) 

k\ \ n / 
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Lemma [52] and (|8.31|) imply 

^ nkz{n) = ^ nkz{n) + ^ nfcz(n) 



Cv2n 



<ne-"z(a) + ^(2BA)('=--')/2 J] z(n') 



and the result for < A; < 2Bn follows by dividing by Z* , with cg = 
1/(25). 

Finally, if > 2Bn, then for every i > 1 we have > i =^ N > 
knk/2 > ki/2 > Bn, and thus by Lemma 18.21 

Hence, Enfc(G;,/„^^) = E£i ^KlG^./^^^) > i) < 2e-^/(2^). □ 

Let Xn^K '■= "-"""^ Z^fcLo ^'^"'fc(^n A/n cs) ^ partial sum of the sum in 
()2.10p for G* x/n-S' '^^6'^) fo'^ every fixed K, by Lemma |8. 




/ k=K+l 



which can be made arbitrarily small by choosing K large. Since X^^k — ^ 
^k=o Po-sIa')!^} as n — > oo for every fixed K, ()2.10p follows by standard 
arguments. (See, for example, the much more general 0, Theorem 4.2].) 
This completes the proof of Theorem 17.31 

9. Proof of Lemma [731 and Theorem 12.11 

Proof of Lemma \7.4\ We use Lemma [7. 21 together with the result of ^] (with 
previous partial results by many authors) that states that, for a sequence of 
degree sequences d = d^") satisfying — > oo, if df = 0(^^ dj) , then 
liminfP(G* is simple) > 0. (The converse holds also, see (§1.) In other 
words, for every K there exist constants ok and bx > such that, if 

(i) E^*-"^ f^^^) E^i-^E^*' (^-^^ 

i i i 

then 

P(Gj; d is simple) > hx- (9.2) 
Let p(d) := P(G* ^ is simple). By Lemma rr2l for every K, 

KG*nM/n;S '^^ ^^ple) = Ep(d(G;,,^/„^5)) 

> hK Hd{Gl^^,^.s) satisfies m) ■ (9-3) 
Thus, it suffices to show that liminf P(d(G* ^^^^.5) satisfies ([9T])) > 0. 
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First, consider the case Jl > 0. By Theorem 17.31 

lT.d^iG:,X^/n■,Sy ^ Ar, T = 1,2, 



for some constants Ar > 0. Hence, taking K := ^2/^11 + 1, 
HdiG:,x^/n;s) satisfies m) - 1 

and the result fohows in this case. 

Now suppose that /2 = 0, which can occur only if G 5. Although the 
graphs are sparser in this case, and intuitively it seems more probable that 
they are simple, we have not found a really simple proof and have to work 
harder in this case. (The proof above is not valid since now Ai = A2 = 0.) 
Let S' := 5\{0}, and let M := n — ng be the number of non-isolated vertices 

Let < m < n and let V be any subset of {1, 2, . . . , n} with \V\ = m. 
If we condition G* . , c on having the set of non-isolated vertices equal to 
V , we evidently get a random multigraph G*^ \„/n-S' °^ ^ (^P *° relabelling 
the vertices) with n — m isolated vertices added. It follows that, for r > 0, 

i i 

Note that the relevant parameter of G* , / is mXn/n, not A„. Since we 
consider < m < n, and the case m = is trivial and thus can be ignored, 
we have < mA„/n < A„ < C15, and thus Lemma 18.71 implies that 

00 

^Y.^^(G*^M/n;S')' =^Y.'''MG*m,X„/n;S') < Gwm, m < U, 

i k=l 

for some constant Ciq not depending on m or n. Furthermore, since ^ S', 
each vertex degree is at least 1 and ^idi{G*^ ^ /n-S'^ — Consequently, 
choosing K = 4Ci6, it follows by Markov's inequality that, for every m < n, 
with probability at least |, 

Y.^^(G*^,x./n;S'f < ^Ciem <KY,d^{Gl^^^^/^.^s')- 

i i 

Consequently, by conditioning on M and using ([97 



MD(ii) holds for > 3 



4- 

Hence, whenever 



then (19. ip fails for G* ^ j^.^ with probability at most ^ + | = |, and thus 

by dasD 

A„/n;5 is simple) > \bK. 
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The only remaining case is when /I = and ()9.5p is false. We recall 
T.idi{G) = Ejjnj{G) and define M* := {n e : Ejj^j < ax}- Thus 
we now have, using (18.20 . 



I < H^(Glx./n;s) e = ^ E ^(")- (9-6) 



^* 

n,\n/n;S neJ\f* 



Further, M* is a finite set (with nj = for j > ax), and (j8.3p yields 
,(„)<4(r„,l,!,f^)'^-'<C„„»-"f^)'^'"', „£A-. (9.7) 



no! V n / v n 

We now use Theorem 14.121 which shows that /I = is possible only when 
1^5. Thus, if n G AA* C Ms, then ni = and i jn^ > = n-no; 

hence, since A„ = 0(1) = 0(n), by (j9.7p . 

/ \ \ n— nn 

z(n)<Ci8n"-"«(^) <Ci9, nGAA*. 



n 



Therefore, by 



neA/"* 



However, if is the empty graph with n vertices and no edges, then by 
», = i?n) = l/^;A„/n;5 > ^2o'- Now is simple, and so 

the result follows in this case too. □ 

Proof of Theorem \2.1\ By Lemmas 17.11 and 17.41 for any event £ and n large 
enough, 

¥{G* . I . G £:) 
'(«-^»/-- ^ ^) ^ P(G:,;*riss..nple) ^ "(«:^-./"^- ^ 

Hence parts (i)-(iii) follow directly from Theorem 17. 31 Similarly, 

^{Gn,\„/n is an 5-graph) = P(G;_;,_^/„ is an 5-graph | is simple) 

^{'^*n\ In is an 5-graph) 
= is simple I is an 5-graph) ^-^^^^^^ ^.^^^^^ 

P(G; ^ , is simple) 

= T^/n* " ' ■ — '■ VT ^(^n,A„/n is an 5-graph) 

^(^n,A„/n IS simple) 



and (iv) follows by Theorem 17.31 and Lemma 17.41 (applied both to S and 



with 5 replaced by Z>o). □ 
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10. Proofs of Theorems 13.11 and 13.61 

Proof of Theorem \3.1[ The case /2 = is trivial by Theorem l2.1l as remarked 
in Remar k 13.21 so we will assume /x > 0. We use the results of Molloy and 
Reed lisjl in the following version, see Janson and Luczak [ll|, Theo- 
rem 2.3 and Remark 2.7]; we only consider the limiting degree distribution 
(Pfc)^o given by pk = Fos{fi){k}. 

Let v = vip), Q = QiJi), 7 and ( be as in Section [3l the existence of a 
unique solution ^ G (0, 1) in (i) follows by [11, Lemma 5.5]. By assumption, 
Po + P2 < 1- Further, let Gn,d be the random graph with given degree 
sequence d, chosen uniformly among all such graphs (assuming that there 

(2) 

is at least one), and let F^^d and F^^ be the largest and second largest 
components of Gn,d- 

Theorem 10.1. Suppose that, for each n, d = (di)^ is a sequence of non- 
negative integers such that Y17=i even, and that 

(i) \{i : di = k}\/n ^ Pk as n ^ oo, for every k > 0; 

(ii) Er=id? = 0(n). 

Then, the following hold for the random graph G„.d; as n ^ oo: 

v(F„,d)/?^ ^ 7: e{Tn,d)/n-^C 

v(rS)/n ^ 0, e(Fi^;j)/n ^ 0. 

This theorem is stated as a limit result, but it can be reformulated as 
follows. 

Theorem 10.2. For every e > and C < oo, there exists 6 > such that 
if n > 6^^ and d = (dj)" is a degree sequence such that Y17=i ^* even and 

(i) Efclo |l{^ '■ di = k}\/n-pk\ < 5, 

(ii) Er=i < Cn, 
then 

IP(|v(F„,d)/ri - 7l > e) < P(|e(F„,d)/ri - CI > e) < e, 

IF'(v(ri'!i)/n >£)<£, IP(e(Fi'i)/n > e) < e. 

By Theorem 12. H for every e > 0, a suitable C and sufficiently large n, 
the random degree sequence d.{Gn Xn/n;s) satisfies the conditions (i) and (ii) 
of Theorem 110.21 with probability at least 1 — e. 

Since (G„,A„/n;5 I ^{Gn,\^/n;S = d) = G„,d by Lemmas O and Ell it 
follows that ¥{\\{Tn,\n/n\s) /''^ ~ l\ > ^) < 2e if n is large enough, and 
similarly for &iXn,\n/n;s) and for F^ ^ /n\S' ^^i^h proves Theorem 13. 1[ □ 

Proof of Theorem 13.61 This proof is similar, using Theorem 2.4]; we 
omit the details while noting that we now need condition (i) of Theorem ll0.21 
and in addition the condition (ii') (stronger than (ii) above) that 'Y^^=i ^ 
Cn for some c > 0. This holds for d(Gn^x^/n;s) with probability > 1 — e for 
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suitable c and C (that may depend on e), as a consequence of the following 
corollary of Lemma 18.71 

Lemma 10.3. Assume that An ^ A > 0. If c < cg, then 

n oo 
i=l k=0 

By Lemmas 17.11 and 17.41 the conclusion of the lemma is valid for Gn^x„/n;S 
also. □ 
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